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We provide a pedagogical introduction to the subject of Standard Model decays of un- 
polarized top quarks into unpolarized and polarized M^-bosons including their QCD and 
electroweak radiative corrections. 

1 Introductory remarks 

These lectures held by one of us (JGK) at the II Hclmholtz International Summer School on 
Heavy Quark Physics in Dubna, Russia (August 11 - 21 2008) are meant as pedagogical lectures 
aimed at the level of the audience which, on the participants' side, was composed of graduate 
students with a few postdoctoral students mixed in. We give many details on the Born level 
calculation of rates and angular decay distributions which can be profitably used in the higher 
order radiative correction calculations. The material collected in the write-up of the lectures 
given by one of us at the International School on Heavy Quark Physics in Dubna, Russia (27 
May - 5 Jun 2002) [1] covering similar topics will not always be repeated. In addition to the 
review [T] we very much recommend the excellent reviews on top quark physics in [51 [3J [31 [51 [S] . 
One of the main aim of these lectures is to illustrate advanced loop techniques in simple Born 
term settings. We begin by listing the basic properties of the Standard Model (SM) top quark 
and its SM decay features. 

1.1 Mass of the top quark 

In our numerical calculations we always take the top quark mass to be rrit = 175 GcV. The 
latest Tevatron combination is TTit — 173.1 ± 0.6(stat.) ± l.l(.s?;si.) GeV [Tj. Since all our results 
are in closed analytical form any other value of the top quark mass can be used as input in 
these formulas. 

There have been suggestions for indirect measurements of the top quark mass through the 
measurements of dynamic quantities that depend on the value of the top quark mass. For 
example, the SM (tf)-production rate at e.g. hadron colliders is sensitive to the value of the top 
quark mass (in particular at the Tevatron II) and thus the (tf)-production rate could be used 
to "measure" the top quark mass. Another possibility is to accurately measure the longitudinal 
and transverse-minus hclicity decay rates of the top quark. The ratio of the two hclicity rates 
is well suited for an indirect determination of the top quark mass since the ratio depends 
quadratically on the top quark mass, i.e. Fl/P^ ~ mj/m^. One should, however, always take 
into account radiative corrections in such indirect top quark mass measurements. For example. 
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in the latter case the NLO QCD and electroweak radiative corrections have different effects on 
the two partial helicity rates which lead to a 3.6% upward shift in the helicity rate ratio Fl/F- 
for top quark masses around 175 GeV [3 [9]. 

In a third method one measures the mean distance that 6-hadrons from (tt)-events travel 
before they decay [10]. The mean distance is obviously correlated with the value of the top 
quark mass. Needless to say that all these indirect top quark mass measurements crucially 
depend on the assumed correctness of the SM. 



1.2 Top quark decays before it can hadronize 

Singly produced top quarks in hadronic collisions are produced by weak interactions and are 
almost 100% polarized. The top quark retains its polarization which it has at birth when it 
decays. The standard argument is that the life time of the top quark (rt ^ 4.6 x IG^^^s) is 
shorter than the hadronization time which is characterized by the inverse of the nonperturbative 
scale of QCD, i.e.Ag^^ = 10~^^s. 

However, one can do better as pointed out in who extended earlier work on depolar- 
ization effects in the bottom sector |121 113) . Consider a polarized top quark which picks up a 
s-wave light antiquark of opposite spin direction. This state will be a coherent superposition of 
the spin and spin 1 mesonic ground states as follows 



The coherent superposition will become decoherent on two counts. First the system oscil- 
lates between the two mass cigenstatcs with a time scale tdocoherencc ~ l/Am^ ~ 6 • 10~^^s 
characterized by the mass difference Am^ = tot* — "m-T ~ {mb/mt)AmB ~ IMeV where 
Ams = niB' ^ niB- Loss of coherence through the decay T* T + j can be neglected since it 
sets in much later at a time scale tdccay ~ 6 • 10~^^s |11| . Thus the depolarization time scale is 
set by tdocohcrcncc and is larger than the tradional estimate based on Ag^^ ^ 10~^'^s by a factor 
of 60. Altogether, the top quark has decayed after tj = 4.6 • lO^^^s much before depolarization 
sets in at tdccohcrcncc ~ 6 • 10~^^s. One concludes that the top quark retains its polarization 
which it has at birth when it decays. 

The decay of polarized top quarks and the corresponding spin-momentum correlations in 
these decays will not be discussed in these lectures. A discussion of the spin-momentum cor- 
relations and their NLO QCD corrections can be found in [S] [T5]. We mention that top 
quarks produced at e'''e~-colliders also possess a high degree of polarization which, in addition, 
can also be attuned by tuning the beam polarization. 

The issue of whether the top quark retains its original polarization when it decays is also 
of importance in the case of hadronically produced top quark pairs. Although the single top 
(or antitop) polarization is zero because parity is conserved in the hadronic production process 
there are sizable spin-spin correlations of the top and antitop quark spins which give important 
information on the (tf)-production process (see e.g. [3]). 




(1) 
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1.3 Dominance of the decay t ^ Xb + 

From the unitarity of the KM-matrix one has the relation 

iKbP + \Vcb\^ + \Vtb\^ = 1 . (2) 

(«0. 004)2 («0.04)2 



which are neghgible compared to the dominant mode t Xb + 



lys 



X, + W+ 



1.4 Rate ratio of t ^ 6 + {-^ leptons) and t ^ h + hadrons) 

Let us list the possible leptonic and hadronic decay modes of the the . For the leptonic 
modes one has the three modes 

W+ [T+Vr), (e+iye) weight : 3 (3) 

When listing the weight factor we have neglected lepton mass effects. 
For the hadronic modes one has 

cb, cs, cd weight : 1 (X) 3 [colour summation) 

ub, us, ud weight : 1 (g) 3 [colour summation) (4) 



using again the unitarity of the KM-matrix X]7=6 s d l^cjP = 1 and X]7=bs(il^iP = 1- 



Again mass effects have been neglected. In ^ we have summed over the respective three modes 

he KM-matrix Y.j=b,s,d l^^^j 
actor of three from colour si 

r(t -> 5 iy+ (-> leptons)) 



■^j—b^s^d I ^Jl iL^j—b^s 

addition one has to add in a factor of three from colour summation. One thus obtains 



T{t^b + W+{-^ hadrons)) 6 



(5) 



1.5 Width of the top quark 



As mentioned before the top quark decays almost 100% to t — > & + W"*" in the SM. The other 
SM decay modes are negligible. Let us list the theoretical values of the SM decay width and 
radiative corrections relative to the Born term width (r(Born) = 1.56 GeV for m?, = 0). 

Tt^^,J^^,+ = 1 Born LO 



r(Born) 

- 0.27% Born mb ^ Q 
-8.5% QCD NLO [16] 

+ 1.55% electroweak NLO [HI [18] 

- 1.56% finite W+ - width [16] 

- 2.25% QCD NNLO [13 [20] (6) 

The NLO and NNLO QCD corrections and the NLO electroweak corections will be discussed 
in Sec. 2. The finite width corrections will be discussed in Sec. 5. 



^In order to simplify the notation we shall in the following refer to the decay £ — > Xf, + iy+ as t — > 6 -|- 
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It is interesting to note that the non-SM decay width into a charged Higgs t — > 6 + can 
become comparable in size to the SM decay width t b+W~^ for small and large values of tan /3 
if 771^+ is not too close to the phase space boundary (see e.g. [4]). A precise measurement of the 
top quark decay width could therefore provide stringent exclusion regions in the (tan/3, m^+)- 
parameter space of Two-Higgs-Doublet models which contain a charged Higgs boson. 

The measurement of the top quark decay width is not simple at hadron colliders. In principle 
there arc two methods to experimentally get a handle on the decay width Tf or lifetime Tf ~ I /Ff 
of the top quark. One can attempt to measure the mean decay length in the laboratory which 
is given by the mean decay length 



Vlab ■ Tlab 



PjC-T = 



Plab 



(7) 



where viab = f3c and Tiab = Jt, and /? = piabc/Eiab and 7 = Eiab / (rnc^) . That this measurement 
is difficult is illustrated by the following example. Take a top quark width of 1.43 GeV. The 
laboratory momentum of the top quark piab must have the astronomically high value of « 
10^^ GeV to produce a mean decay length of 1mm. Nevertheless CDF has attempted such 
a measurement using information on the magnitude of the impact parameter of the charged 
lepton with respect to the collision vertex. CDF puts a 95% confidence level upper limit of 
1.8 • 10~^^s on the lifetime of the top quark which corresponds to a 3.7 • 10~^^ GeV lower limit 
on the top quark width [21]. Naturally, this is not a very useful bound. CDF also provides an 
upper limit on the top quark width by a fit of the reconstructed top quark mass to a Breit- 
Wigner shape function. The reult is Ft < 13.1 GeV at 95% C.L. [22]. The upper bound is still 
nine times larger than the expected SM width. 

An indirect way of determining the top quark width relies heavily on the validity of the 
SM. The suggestion is to measure the branching ratio B{t bW) = F(i bW)/r{t all). 
This could be done e.g. by measuring the rate of top quark pair production followed by their 
decays t bW, i.e. by measuring (Tu ■ {B{t bW)Y . Assuming that one can reliably 
calculate ati one can then extract B{t — *■ bW) (see e.g. |3]). In the simplest version of this 
approach one takes the SM value V{t ^ bW) to determine the width of the top quark through 
T{t all) = T{t bW)/B{t bW). In a more sophisticated approach one uses single-top 
production to extract the parameters that determine the partial width F(i bW) [23] . 

We mention that a much improved determination of the top quark width with an uncertainty 
of AF w 30 MeV can be expected from a multi-parameter scan of the threshold region of (tt)- 
production at the ILC [M] . 



1.6 Top quark yield 

At the LHC top quark pairs will be produced quite copiously in 7 on 7 TeV proton-proton 
collisions. After a one-year probation run at reduced energies and luminosities starting in the 
end of 2009 the LHC will start running at full energy in 2010 with a low luminosity run of 
L « 10^^cm~^s~"'^. After a luminosity upgrade around the year 2017 the high luminosity run 
will have L w 10^'*cto~^s~^. Multiply these numbers with a{tt) ^ 825 p6 = 825 x 10~^^cm^ to 
obtain w 1 (10) (tf)-pairs every second for the low (high) luminosity run. 

Top quark pair production at the Tevatron II (1 on 1 TeV pp-collisions) occurs at a reduced 
rate. Because the energy of the Tevatron II is lower, the (tf)-production cross section is down 

■^We have employed a mixed notation in the last equality of Eq. Jt} where we set c = 1 for the quantities in 
the square bracket [piab/n^]. 
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by a factor of ss 100. In addition, the Tevatron II luminosity is down by a factor of « 10 
compared to the LHC low luminosity run. Taking these two factors into account one has 
-pairs per second at the Tevatron II. 

In the SM single top production cross section in (pp)- and (pp)-collisions is down by a factor 
of w 3 compared to top quark pair production. The nice feature of single top production is 
that the top quarks are polarized since the production of single top quarks proceeds through 
weak interactions. The polarization can be calculated to be close to 100% (see e.g. [4]). 

At the ILC (tF)-production will occur at a somewhat reduced rate compared to the LHC. 
At 500 GeV the NLO rate is a{tt) - 0.5 pfe = 0.5 x IQ-^'^cm'^ (see e.g. [23 ES]) which gives 
10~^ (tt)-events per second assuming a luminosity of i w 2 • lO'^'^cm^^s"-'^. 



1.7 Polarization of gauge boson 

The decay t ~> b + is weak and therefore the VF+-boson is in general expected to be 
polarized. We shall refer to the three partial rates that correspond to the three polarization 
states of the W^+-boson as longitudinal (F^), transverse-plus (r+) and transverse- minus (r_). 
At leading order (LO) the results for the helicity fractions Gi — r^/F (i = ) (or, in 

another language, for the normalized diagonal density matrix elements of the M^^-boson poO) 
and p ) ar^ 

1 2y2 

where = m'^/m^ = 0.211 with mt ~ 0. Numerically one has 

Gl-G+'-G-^ 0.703 : : 0.297 . (9) 
Note that Gl + G+ + G- = 1. In comparison, an impolarized would correspond to 

Gl:G+:G- = 1/3 : 1/3 : 1/3 . (10) 



1.8 Dominance of the longitudinal mode 

As (myy/'m-t) the longitudinal polarization vector becomes increasingly parallel to (see 
e.g. [U), viz. 

e1^—(q'' + 0{mw/mt)). (11) 

Therefore the longitudinal mode dominates in the large top quark mass limit. In fact, from 
g^Ufc7^(l — 75)ut = 77itUfc(l-|-75)ut one concludes from dimensional arguments that ^ Gpra\ 
whereas F± ~ GpiTitm^. or T±/Tl ~ m^/m^. 

An explicit calculation shows that F_|_ = at LO for mf, = (see Eq.®). Looking at Fig. 1 
the vanishing of the LO transverse-plus rate F_|_ can be understood from angular momentum 
conservation. First remember that a massless left-chiral fermion is left-handed as drawn in 
Fig. 1. At LO one has a back-to-back decay configuration. Therefore the iy+-boson cannot 
be right-handed because the TO-quantum numbers in the final state would add up to 3/2 which 
cannot be reached by the spin 1/2 top quark in the initial state. At NLO (or any higher order) 
the decay t ^ b + g + is, in general, no longer back-to-back as illustrated in Fig. 1 and one 

^The helicities of the W-boson are alternatively labelled by (L, — ), (0, -1-1, —1) or by (L, T+ , T_ ). 
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anticipates that r+ ^ at NLO and at any higher order. This is, in fact borne out by the NLO 
calculation to be described later on. The physics interest lies in the fact that nonvanishing 
transverse-plus helicity rates can also be generated by non-SM right-chiral (t&)-currents. In 
order to unambigously identify non-SM contributions to the transverse-plus helicity rate it is 
therefore important to get a quantitative handle on the size of the SM higher order radiative 
correction contributions to the transverse-plus helicity rate. 




Figure 1: Angular momentum conservation for t b + and for t b + W'^ + g. 



1.9 Measurement of the helicity fractions of the through the 
angular decay distribution in its decay 

The decays weakly to {l^vi) or to (qiqj)- The angular decay distribution can therefore be 
utilized to analyze the polarization of the decaying W~^, i.e. the W'^ is self-analyzing. 

The has the three (diagonal) polarization states L, and r_ the weights of which are 
determined by the three partial helicity rates Tl and r± . As we shall explicitly derive further 
on, the angular decay distribution for t ^ b + W^{-^ + I'l) reads 




Figure 2: Definition of polar angle 6 in the rest system. 

where the polar angle 6 is measured in the M^-rest frame as shown in Fig. 2. Integrating over 
cos 6 one recovers the total rate F = F^ -|- F+ -|- F_ . If the were unpolarized one would 
have Fi = F+ = r_ = F/3 resulting in a flat decay distribution dT/dcos9 = F/2 . 
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One can also define a forward-backward asymmetry by considering the rate in the forward 
hemisphere Tp and in the backward hemisphere in the M^~''-rest frame. The forward- 
backward asymmetry Afb is then given by 



Afb 



Tf-Tb 



3 r+ -r_ 

4 Fl + r+ + r 



(13) 



At the Born term level one has 



Afb (Born) 



3 2y^ 

4 l + 2y2 



= -0.22 . 



(14) 



The forward-backward asymmetry is negative, i.e. one has more leptons in the backward 
hemisphere than in the forward hemisphere. The numerical value of the forward-backward 
asymmetry is not very large on account of the dominance of the longitudinal mode. 

It is always useful to check on the correctness of the sign of the parity violating term pro- 
portional to (±cos0) and thereby on the sign of ApB- This is again easily done by considering 
the collincar cases cos9 — ±1 and appealing to angular momentum conservation. And, in fact, 
Eg. lfT^ shows that the mode r_ decouples in the forward direction cos 6 — 1 (and vice versa r_|_ 
decouples in the backward direction) as can be appreciated from the helicity configurations in 
Fig. 2. This implies that r+ favours forward leptons leading to energetic leptons in the i-rest 
frame whereas r_ favours backward leading to less energetic leptons in the t-rcst frame. As 
we have seen r+ = at LO so that one expects a softer lepton spectrum in the t-rest frame 
then in the case of the decay of an unpolarizcd . 

2 Top quark decay rate 
2.1 Leading (LO) rate 

We shall calculate the leading order rate in three different ways for pedagogical reasons. The 
first way is the traditional covariant way where no particular sophistication is needed. In the 
second way we use the helicity methods which has the advantage that by calculating the helicity 
amplitudes one has the full spin information of the decay at hand. In the third method we use 
the optical theorem which serves the purpose of introducing rather sophisticated technical 
material in a simple setting which are needed later on in the higher order calculations. 

2.1.1 Covariant method 

The matrix element for the decay t b + {pt = Pt + l) is given by 



M = 



■ 9w 



(15) 



2V2 



Upon squaring and summing over the spins one obtains 



(16) 



spins 
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where we write Jospins l-^l^ ~ Use of the completeness relations 

= |1 + TO (17) 



±1/2 



and 



leads to {mi, — 0) 



^e''(m)e-(TO) = -g^'' + ^ (18) 



0,± 



2 

|M|2 =|i|14,|2Tr {|^^7^(1 - 75)(^, + m07^(l " 75)} ( - 



--%Vtb\H({ptPb) + -^2{p,pt) {ptq)) . (19) 



Using four-momentum conservation pt = Pb + Q and the mass shell conditions p^ = = 
and = to^ one obtains 



|Mp = 4l^*&l'4TO?l^(l + 2y2) 

y 



The rate can be computed using the two-body decay formula 

where R2 denotes the two-body phase space integral ^7\. We symbolically write R2 [|M for 
the two-body phase space integration over the squared matrix clement |Af i.e. we write 

[ ^ i / (5^? 2^ / (2^?^ p"' - - i»r (21) 

In order to stay general we calculate R2 for mi, ^ 0. The phase space integral will be evaluated 
in the top quark rest system. We write l/{2Ew) = J dEw 5{q^ — m^) = J dEw 5{E'y^ — 
— TO^), where we implicitly take the positive energy solution Ew = +\/|g|^ + The 
corresponding relation for the bottom quark energy reads l/(2Ei,) = J dEi,d{pl — to^). Using 
these two relations one converts the three-dimensional integrations in (j2ip into four-dimcnsional 
integrations. One obtains 

R2 [ \M\'] = J d\ J d% S{q' - ml,) 5{pI - ml) (p, ^ p, - q) |M|2 . (22) 

The integration over d^pb can be done with the result that the argument of the second (5-function 
becomes [pt — q)^ — ml = TOj — 2mtEw + to^ — to^, i.e. 

R2 [ \M\^] = j dS S{q' - m^^) 6{ipt- q)' - ml) \Mf . (23) 
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Next one integrates over dEw = d{2mtEw) /2mt with the result that the argument of the 
remaining 5- function becomes — = E'^ — \(}\^ — ~^ ("^t + "^VF ~ "^6)^/(4^4) — 
[(fl^ — iTi^. The remaining integration over d'^q can be done using spherical coordinates such 
that d^q dn\q\^d\q\^ = ^\Mq\^. The result is 

R^[m']=^^\<l\\M\\ (24) 

where \q\ — \/ X{mf, m^, rn^) / {2mt) is the magnitude of the momentum of the M^+-boson in 
the top quark rest (^-rest) frame and where A(a, 6, c) = (a^ + fe^ + — 2a6 — 2ac— 26c) is Kallen's 
function. Naturally we could have calculated the two-body phase space R2 directly without 
including the squared matrix element |Mp in the integrand as long as the kinematic variables 
in |Mp arc fixed according to four-momentum conservation and the mass-shell conditions. 

We now return to the approximation nri}, = where \q \ — mt{\ — ij^) /2. Substituting the 
matrix element squared (fTQ]) into the rate formula ([20|) one obtains 



r(Born) = To (1-2/^)2(1 + 22/2), (25) 
where Fq is the mw = Born term rate [g^ / (8m^ ) = /^/2) 

2.1.2 Helicity amplitude method 

The helicity amplitudes for i ^ 6 -I- VF^ can be calculated from the transition matrix element 
by using spinors and polarization vectors with definite helicities At, A{, and \w- One needs to 
calculate (we omit the coupling factor ~*^^^*b) 

Hx,..x,x,, = ub{\)r{l - ^5)ut{Xt)e;{Xw) ■ (27) 

We shall work in the i-rest system with the z-axis along the (see Fig. 3) such that At = 
—Ah -I- Xw- In order to be general we keep mi, ^ 0. 



-AAAAAAAAy W 
t 



Figure 3: Definition of the the two-body coordinate system in the top quark rest system. 



Let us collect the relevant t-rest system spinor and polarization vector expressions. For the 
helicity spinors one has 



itt(l/2) =V2^ 
itt(-l/2) =V2^t 



(^q) , "^(1/2) = /B^T^ ( E^X- ) 



X- 




(-1/2) = ^/WT 



■nib 



(28) 
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where x± are Pauli spinors given by x+ = (q) SiXid x- = (i) 
The hehcity polarization four-vectors of the W'^ read 



1 



e*(±l) = — (0;±l,-*,0) 



4(0) 



V2 
1 



(|(7l;0,0,-Qo) . 



(29) 



There are altogether four possible helicity configurations in t ^ b + which are listed in 
Table 1. 



At 




Xw 


1/2 


-1/2 





-1/2 


1/2 





1/2 


1/2 


1 


-1/2 


-1/2 


-1 



Table 1: Helicity configurations in t ^ b + W^. 
For the helicity amplitudes H\^.x^\„- {Q± ~ {rrit ± m^)^ — q'^) one obtains 



Hi.i 

2 ' 2 



H_ 



= -mt{^ 

1 = -V2(. 
1 = -V2(. 



{+ - 

?7 



TTlf,— »-0 



- 2m2v/l-2;2 

0, 
0, 



-2V2mt^/l'y\ (30) 



where we have included both the nif, ^ and mf, = results in (j30p . The squared matrix 
element |Mp finally is given by 



At — — Afa + Aw 



4m: 



2(1 -J/') 






T+ 



(31) 



where we have set = in the second line of ([?!]) . The result agrees with the covariant calcu- 
lation (see Eg. ([TO]) ). The advantage of the helicity method is that one can separately identify the 
three (diagonal) helicity contributions of the boson L, r+ and T_ as indicated in Eg. pTj) . 
In fact, the helicity amplitudes contain the complete spin information of the process. Thus one 
can easily calculate other polarization effects using the helicity amplitudes such as the decay of 
polarized top quarks, the polarization of the bottom quark and polarization correlation effects, 
mb ^ effects arc easily included by using the mfj ^ helicity amplitudes in Eq. (|30p . One can 
also define covariant helicity projectors which allow one to directly calculate the longitudinal, 
the transverse-plus and transverse-minus helicity rates without taking recourse to the helicity 
amplitudes. This will be described in Sec. 5. 
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2.1.3 Optical theorem method and cutting rules 

In this subsection we shaU use yet another method to calculate the leading order rate for 
t ^ b + using the optical theorem. Whereas the optical theorem method does not offer 
particular technical advantages in LO calculations it is the method of choice for higher order 
calculations as e.g. the calculation of the NNLO rate to be described later on. The reason is 
simply that the phase space integrations in the NNLO radiative correction calculations become 
prohibitively complicated and cannot be automated as easily as higher order loop calculations. 
We present the optical theorem method for the LO case for pedagogical reasons because the 
LO discussion allows us to introduce concepts which are also needed in the NNLO radiative 
correction calculation to be described later on. 




Figure 4: Illustration of the optical theorem method to calculate the LO top quark width. 



The optical theorem relates the width F of a particle to the imaginary part of the self-energy 
contribution S of the particle. In the top quark case one has|f| 

1 ImS 

2st + l mt ^ ' 

where, for the present purposes, S is the one-loop self-energy of the top quark as illustrated in 
Fig.4. 

Using standard Feynman rules [27l [28] the one-loop self-energy contribution is given by 

{pt =Pb-q) 



E f d a I n... 1 — \ «: 

u{Pt,St) J 



i—7=T — — Vtb I o 



(27r)4 VV2 ' 2 pI 

One can again use the completeness relation X]±i/2 ^ + '^0 rewrite Eq. ([33|) as a 
trace. The trace can be taken as in Eq.p9|) except that one now cannot avail of the mass-shell 
conditions = and = 0. One obtains 

V „l5^|T. ,2 8 /• d^q ptq{2ptq + 2q^ + ml,) + mjml, 

^one-loop — --Z-\Vtb\ — — / TTT^ "jV 1 \7 , -IF? 2 , • 1 ■ V"J^J 



J (27^)^ [{pt + qY +it][q' -ml, + ie] 

The usual procedure is to expand the g-dependent numerator factors in terms of the g-dependent 

= cY - m^, 



denominator factors Dq = q^ — m^, and Di, — (pt + g)^ in order to obtain g-indepcndent 



very nice discussion of the optical theorem and related technical material relevant to top quark decays 
can be found in the thesis of I.R. Blokland 1281. 
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numerator factors (the corresponding integrals are called scalar integrals) after cancellation. 
We therefore write 

Ptq = -^{m'f + m'^ + Dq - Db) , 

D, . (35) 



^2 _ 2 



m 



w 

The contributions proportional to Dq and Db cancel against the denominator pole factors and 
their contributions can be dropped when taking the imaginary part since single or zero pole 
contributions have no imaginary part (see e.g. [28j). We therefore have 

^one-loo, = \'^\V,b? S^(l - V'){1 + 2y') I , ^. , -1?. ^ • (36) 



I ml,' ' ' ' J (2^)4 [(p, + q)2 + ie] [g2 _ + 

According to the cutting rules the discontinuity of a Feynman graph is obtained by the 
product of the discontinuities of the pole factors which are being cut, where the discontinuity 
of a single pole is given by [27l [28] 



Disc — \ = -2m S(p^ - m^) . (37) 

Furthermore, the imaginary part and the discontinuity of a graph M are related by 2i ImAf = 
DiscM. One therefore has 

ImS__,„„, = il|i|l^,,|2^(lV)(l+2y') 1 -^S-^mfb{(p,^qf)8{q^-m\,). (38) 

In order to exhibit the similarity to the integral we change the integration variable q — * —q. 
One can then use the result of Sec. 2.1.1 

d\8{q^ ~ ml,) 5{{pt - qf) = |(1 - y^)) 
to arrive at 



where, as before. 



r(Born) =ro(l-y2)2(i + 2y2), (39) 

Fo = -—j=-\Vtb\ and y = . 40) 

8V27r mt 

As it must be the result agrees with the covariant and helicity amplitude calculations. It is 
quite reassuring that the decay rate turns out to be positive definite in the end, as it must be, 
considering all the minus signs and the factors of (i) appearing in the rate calculation using the 
optical theorem method. 



2.1.4 Expansion by regions and the (mvi//'Tit)-expansion 

In Sec 2.1.3 we have calculated the leading order rate by using the optical theorem and cutting 
rules to determine the imaginary part of the one-loop self energy diagram. In this subsection we 
shall go one step further and calculate the leading order rate using a {mw /mt)-expansion which 
allows us to introduce the concepts of expansion by regions and integration-by-parts identities. 
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All latter three concepts are essential in the calculation of the NLLO rate presented in [191 120] • 
As emphasized before we shall pattern the LO rate calculation after the NNLO calculation 
entirely for pedagogical reasons. In the LO case the follow-up calculations are simple enough to 
be presented in a few simple lines, whereas they are more involved in the ftiU NNLO calculation. 

Let us summarize the main ideas of the NNLO rate calculation presented in [l9l [20] which 
we down-size to the present LO case. 

• Reduce the two- mass-scale problem {mt,mw) to a one-mass-scale problem (mt) by ex- 
panding in the ratio mw/mt- Obtain the results as an expansion in powers of mw/mt. 

• Use dimensional regularization to regularize the UV and IR/M singularities 

• Use the method of expansion by regions to calculate the one- loop integral [29l [30l [31] . 
One has to consider the two regions [29l [30l [31] : 

• Hard region 

The loop momentum is hard and is of 0(rnt). One can then expand the VF-propagator 
as a power series in m^^/rn^ ^ 1: 

^-L^=i+4^+4^+...=if:f4^V • (41) 

g2 - (?2 ^4 ^6 q2 Z^^y q2 J 

The massive propagator has thereby been converted into a sum of massless propagators. 

• Soft region 

The momentum q flowing through the W is soft. One therefore cannot use the above 
expansion ([1T|) of the W propagator. However, in the soft region one can expand the 
6-quark propagator, cif. 

{Pt + q? Pits P" ) ■ 

There is only one denominator factor in the loop integral and its imaginary part vanishes 

Therefore there is no contribution from the soft region in the one-loop case. This is 
different at NLO and NNLO. 

What remains to be done is to evaluate integrals of the form 

Inji^ I (44) 

which result from the / expansion in the hard region. The integrals can all be reduced 
to one master integral by using integration-by-parts identities. 

The first term in the expansion (|4ip . n = 0, leads to a two-point one-loop integral of the 
form 

1 d^q l_ 

I (2^)4 q^{pt + qY 



/ -tt^^.-tt:—^, ■ (45) 
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We calculate the one-loop integral directly in dimensional regularization {D = A — It) and 
take its imaginary part at the end without resorting to the cutting rules. The details of how 
to evaluate one-loop integrals in dimensional regularization can be found in [27] ■ One first 
introduces a one parameter Feynman parametrization, collects terms and performs a shift in 
the integration variable (q -I- xpt) — > g, i.e. 



1 



dx- 



{2^)D q^{pt + qY J {2n)D [{q + pt^ x + {1 - x)f 

d^q , 1 



dx- 



(27r)^yo [{q + xptY + pM'i^ - x)Y 

t dx ^ 2 . (46) 

(27r)^io [q' + pix{l ~ x)f 

Next we do a Wick rotation go ^ iqoE- The factor of i from the Wick rotation cancels the 
factor of i in the denominator of ()44p . One then does a D-dimensional Euclidean integration 
over the loop momentum g, and, finally, one integrates over the Feynman parameter x which 
results in Euler's Beta function B{\ — e, 1 — e). The sequence of steps is represented in the 
following sequence of equations: 

Im / — - — TT^—^, = Im / dx- 2 ) 



i {2tt)D q-^{pt + qY 7o (47r)^/2 r(2) \-plx{l - x 



1 



(4^)D/2 r(l + e)^ 

r(l + 6) T{e) T\l-e) 
(4^)D/2 r(i + e) r(2 - 2e) 

We retain only the finite term in the last line of (|T7l) . One obtains 



(mj)"' siuTre . (47) 



We have used 



which leads to 



1 d-^g 1 _ 1 

^™ ' 7 (27r)4 q^{p.t + g)2 = 16^ ■ ("^^^ 



2 

^ ein(-p?/m?)-^ ^ e-^in(-p?/m?) ^ - eln(pV"^?) + •■• (49) 



lm(^)-' = lm(-l + iO)-' = sinTre = Tre-f ... . (50) 



m 



H 

In addition to the integral (|45p with n = the imaginary part of which we have just 
calculated wc also need the imaginary parts of the integrals (|44|l with n > 1. They can be 
obtained from the "master integral" (|45p by integration- by-parts (IBP) techniques [521 [33]. The 
general procedure of reducing a set of integrals to a set of simpler integrals is called "reduction 
to master integrals". In the present case this reduction is quite trivial but can become quite 
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involved in more general settings. The reduction procedure has been automated by the Laporta 
algorithm [Ml [35]. 

Technical aside: Integration-by-parts (IBP) identities [5^ [55] . 

In order to calculate the integral corresponding to the second term in the expansion (|4ip we 
consider the differential form (5^ := d/dq'^) 

^ [pt + qY _ d^ipt + qY , [pt + qY ^ 1 ^ {Pt + qY ^ 1 ... x 



""q^iPt + q? q'{Pt+qY q^ {Pt + qY {Pt + q? ^ q 

Differentiate carefully, i.e. d^j^q'^ :~ = D = 4 — 2e, and drop the "surface term" on the 
left-hand side. Also use 2ptq ~ ^itl^ — q^ + {pt + qY- This gives 



1 _ 1 / 2e - 1 1 
qHjhTqY ^ m! Wipt + qY ^ q^ 



(52) 



In dimensional regularization massless tadpole (single pole) diagrams are zero, i.e. one can drop 
the second term on the r.h.s. of ([5^ after dimensional integration. At the relevant order of e 
one therefore has 



I i (47r)^ q^{pt + qY m^t ^"'/ ^ (4^)^ qHPt + qY ' ^^^^ 



Going through the same exercise for ^^'iy^lt^J^^^^yi for n> 2 one finds 

f 1 d^a 1 

/ -JT^ , 4-1/ ^ = . for n>2 (54) 

Because the higher order terms vanish we only need to sum the first two terms in the 
expansion (|4ip . The result 



Im 



i (47r)^ ipt + qY q^ q^ J 167r2 



is in agreement with the one in Sec. 2.1.3. 



2.2 Next-to-leading order (NLO) QCD corrections 

The traditional technique used for NLO calculation is to calculate the one-loop and tree-graph 
contributions separately. In the present case the UV singularities are regularized by dimen- 
sional regularization whereas the IR/M singularities are regularized by introducing giuon and 
bottom quark masses. The IR/M singularities will eventually appear as (In mg)— and (Inm;,) — 
singularities and cancel among the one-loop and tree graph contributions [3 [9] [M] [TH] [36] . We 
mention that the calculation can also be done in dimensional regularization without recourse 
to the traditional Trig ^ and ^ regularization [37j . 

For example, generic diagrams for the QCD NLO calculation are displayed in Fig. 5. 
Without going into the details of the calculation (see e.g. (5) we just quote the result of the 
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1-loop ® Born 



tree (g) tree 





Figure 5: Generic NLO QCD contributions 



NLO calculation. For the total rate one obtains (f = r/r(Born)) (see e.g. [36) ) 

t(Nrn\ = / (l-y^)(5 + V-6y^) 2il-yyil + 2y')n^ 

^ ^> +27r ^(l-y2)2(i + 2y2) I 2y2 3y2 



(1-^2)2(5 + 4^2) 4(l-y2)2( l + 2y2) 

,2 ^ ,,2 



ln(l-y2)_Ii^J^^4iZ±^ln(y)ln(l-y2)_4(l+y2)x 



y J 

The numerical value of the NLO QCD correction appears in Eq.®. Our numerical input 
values are nit = 175 GeV and mw = 80.419 GeV. The strong couphng constant has been 
evolved from as{Mz) = 0.1175 to as(mt) = 0.1070 using two-loop running. Numerically one 
has r = r(Born)(l — 8.54%). One sees that the NLO QCD corrections reduce the Born term 
rate by the large amount of 8.5%. 

In the limit y — *■ one obtains 

f(NLO)^l + ^C.{|-H,^},. (57) 

The leading y — * contribution reduces the rate by 9.26% which is already quite close to the 
rate reduction of the full result (8.5%). We shall return to an assessment of the quality of the 
2/-expansion later on. It is curious to note that the radiative QCD corrections reduce the LO 
rate whereas the radiative QCD corrections to the decay Z ^ qq enhance the LO rate ratio 
f (LO) by fj(6/4)CF = a^Tr, i.e. f(NLO; Z ^ qq) = 1 + as/n. 

The NLO rate can also be calculated by the optical theorem method using the y-expansion. 
At NLO one has contributions both from the soft and the hard region leading to an infinite 
power series in y and ylny where the (y In j/)-contributions come from the interplay of the soft 
and hard integration regions. The results of the y-expansion have been checked against the 
exact result Eq.^ up to 0(y^^) [3H] (see also [25]). 



2.3 NLO electroweak corrections 

In Fig. 6 we have drawn the LO diagram and the four NLO tree-level diagrams that contribute 
to t ^ b + + (7). We use the Feynman-'tHooft gauge so that one has a NLO contribution 
from the charged unphysical Higgs boson as shown in Fig. 6. Compare the number of 
four electroweak NLO tree-level diagrams with the two QCD NLO tree-level diagrams. When 
squaring the tree-level diagrams one would expect a four-fold complexity factor when going 
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Figure 6: Born and electroweak tree-graph contributions to t b + (7). denotes the 
charged Goldstone boson. 



from QCD to the electroweak tree-graph corrections. It is therefore quite remarkable that the 
squared tree graph expressions in both cases are similar in length and structure [S] ■ 



Figure 7: Eighteen electroweak three-point one-loop graphs in Feynman-'tHooft gauge con- 
tributing tot^b + W+. 

In addition to the tree graph contributions one has to consider 18 three-point one-loop 
graphs in the Feynman-'tHooft gauge as shown in Fig. 7. Looking at Fig. 7 one would superfi- 
cially expect 4-1-8-1-8=20 one-loop contributions. However, since there is no (Vl^"'" W^+x°)-vertex, 
this number reduces to 18 as stated before. In Fig. 7 and x° are the charged and neutral 
unphysical Goldstone bosons, and H is the physical Higgs. The results of calculating the one- 
loop contributions exist in amplitude form [17]. In the course of calculating the electroweak 
radiative corrections to the partial helicity rates the results of [T7] were recalculated and con- 
firmed by us. In particular we checked the results of [T7] numerically with the automated loop 
calculation program XLOOPS/GiNaC developed at the University of Mainz [321 HOI [H] . In 
addition to the one-loop three-point functions one has a large number of one-loop two-point 
functions needed in the one-loop renormalization program. Again these have been reevaluated 
using XLOOPS/GiNaC. 

We have used the so-called GF^renormalization scheme for the electroweak corrections where 
Gf, Mw and Mz are used as input parameters. The Gj^-scheme is the appropiate renormal- 
ization scheme for processes with mass scales that are much larger than Mw as in the present 
case. The electroweak radiative corrections are substantially larger in the so-called a-scheme 
where a, Gp and Mz are used as input parameters. The numerical results of the electroweak 
corrections to the rate are given in Eq. ([6]) . 



In the NNLO case squaring of the contributing tree and loop diagrams leads to the four generic 
contributions shown in Fig. 8. However, with present techniques, this method is not viable, 




2.4 NNLO QCD corrections 
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2-loop (8 Born 



1-loop (g) 1-loop 





1-loop gluon emission (g) 1 gluon emission tree (g) tree 

Figure 8: Generic NNLO QCD contributions. 



mainly because the NNLO phase space integration become too difhcult. 

Instead, one resorts again to the optical theorem and calculates the NNLO rate from the 
three-loop self-energy diagrams according to [TH 



r(NNLO) = 



1 



1 



Iml](3 — loop) , 



(58) 



2st + Irrit 

There are altogether 38 three-loop Feynman diagrams a sample of which are shown in Fig. 9 



(a) 



id) 




Figure 9: Sample three- loop diagrams whose imaginary parts contribute to the NNLO calcula- 
tion of the top quark width. 



The main ideas of the NNLO calculation of the rate have already been described in the 
calculation of the Born term rate in Sec. 2.1.3. It turns out that again one only has to consider 
two momentum regions. In the hard region all loop momenta are hard and the VK-propagator 
can be expanded into a series of massless propagators as in the LO case. In the soft region 
the gluon momenta are hard but the loop momentum flowing through the W is soft. Differing 
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from the LO calculation one now also has contributions from the soft region. In the soft region 
the integrals factorize into two-loop self-energy-type integrals and a one-loop vacuum bubble 
diagram which are not difficult to integrate. The interplay of the hard and the soft region leads 
to additional (y" In j/)-terms in the y = (77ivv/mt)-expansion. 

One can reduce all integrals to 23 master integrals by integration- by-parts identities. Use 
was made of Laporta's algorithm in this reduction to master integrals. The imaginary parts of 
the master integrals were calculated using the cutting rules where care had to be taken that 
some of the master integrals admitted several ways of cutting them. We mention that the 
calculation had been done in the general covariant gauge —g^'^ + — \)k^^k'^ /k'^ for the gluon 
in order to check on gauge invariance. The numerical results on the NNLO QCD corrections 
are given in Eq.®. 

3 Vi/-helicity fractions in top quark decays 

3.1 Angular decay distribution for t ^ h + > -|- z/^) (J) 

In Fig. 6 we display the LO amplitude contribution io t ^ b -\- + vi. On squaring the 
amplitude and taking the spin sums one is led to the contraction L^^H^'^ (JioTm). For the 




b 



Figure 10: LO Born term contribution to i ^ & + W^{-^ t'^ + vg) 
lepton tensor we obtain 

L^''= ^Tr|^^7^(l-75V,7''(l-75) 

= P^,Pl+p'iP^u-\^W -^^'""'^'p^o.Pup . (59) 
The LO hadron tensor is given by (m^ = 0) 

H^^^iBoin) = ^TT{f^ + mt)7^(l - 75V,7''(1 " 75) 

= p^pl + p^.pI - pt ■ p.g^'' - ie^'"^'^'p,^,p,0, . (60) 

The factors 1/8 have been introduced for convenience. The result of contracting the lepton and 
hadron tensor reads 

L^,H^''(Born) = 4(pt • pi) {pt ■ p.) = 4(pt • pi) {pt -{q-pe)) ■ (61) 

Note that one originally had 
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L,,.iJ^'^'(Born)(-g;:, + + ^) (62) 



which turns into L^^H^^'^ (Qoixi) in the zero lepton mass case where q^L^'^ = Qi^L^"^ = 0. The 
lepton mass corrections arc of 0{m?^/raf) and arc thus ncgUgiblc. If one wants to include lepton 
mass effects one has to retain the full Vt^-projector in ((62|) . 



One must evaluate the invariant L^^H^'^ in one frame. Here we choose the rest frame of 
the top quark. Since we want to evaluate L^yH^^^ in terms of the angle cos0 defined in the 
W^+-rest frame {Wr.f.) as shown in Fig. 2 wc writ^ 

P£ (W^r.f.) = ^^^^(I;sin6'cos(/),sin6'sin0,cos6') . (63) 

We then boost the lepton momentum (Wi-.f.) to the top quark rest frame (tr.f.) where the 
invariants in (1611) are to be evaluated. The relevant Lorentz boost matrix reads 



L{boost) = n n _ n (64) 





( 90 








k1 \ 


1 





mw 
















raw 







^ k1 








90 / 



such that 

(tr.f.) = i(boOSi) pI (Wr.f.) . 

The boost will not affect the transverse components /x = 1, 2 but only the zero and longitudinal 
components ^ = 0, 3. In Ea. (j64p and |g| denote the energy and momentum of the VF-boson 
in the top quark rest frame. 

In the following we set mi, ~ such that go = ^(1 + J/^) a-^^d \q\ = ^(1 — v'^)- Boosting 
Pi (Wr.f.) one obtains 

p^l{t,,i,) = ^((l + y2) + (l-2/2)cos6'; 2y sinflcosc/), 2ysin6'sin0, (l-2/2) + (l + y2-) ^^g^^ 

The remaining momentum four-vectors in the t-rest frame are given by 

pt = mt(l; 0,0,0), 

V', = ^(l-y^;0,0,-(l-y^)), 

= ^(l+2;^0,0,l-y2)_ (gg) 

We are now in the position to evaluate the invariants appearing in Eq. (|6ip . We sort the 
resulting expression in terms of the polar angle factors sin^ 9 and (1 ± cos6')^/2. Since we are 



^ In Eq.| |63I | we have specified the azimuthal dependence of (VFr f ). This is not really needed in the present 
application because we do not specify a preferred transverse direction. In general, a transverse direction could 
be defined by the polarization of the top quark or the decay products of the 6-quark. In this case one has to 
retain the azimuthal dependence of the lepton's momentum as done in l|63| l. Whereas the sign of polar angle 
correlations can always be checked by physics arguments, there arc no ready physics arguments to check the 
signs of the azimuthal correlations. To get the signs of the azimuthal correlations right it is indispensable to use 
the boosting method as described above (see e.g. 1421 ). 
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not interested in the azimuthal angle dependence in the present application we integrate over 
the azimuthal angle 0. One then obtains the angular decay distribution 



2"" Q ^4 




i,.H^''(Born) = 2.^^| \{1 - y^) \sin^ 9 

^ ^(1+COS0)2 

^(1- COS 0)2 |, (67) 

where, by comparison with Eq. (j2ip . we have identified the three LO hadron contributions 
proportional to L and T±. The normalized helicity fractions Ql and Q± written down before 
in Eq.® can be read off from Eq. (|67p . As we shall sec later on from an angular momentum 
analysis, the sorting of the angular contributions in (|67|) should be done exactly along the three 
angular factors proportional to sit? 9 and (1 ± cos^^)2/2 discussed above. The corresponding 
coefRcient factors arc then proportional to the partial helicity rates and F-t, respectively. 
An untreated and unsortcd Mathcmatica output of Lf^,^H^^ would, in general, lead to quite 
lengthy and messy expressions. 

Repeating the same exercise for mb ^ one obtains {x = mb/nit) 

L,.i/'-(Born) = 27r?^( Uil - x^ - y^l + x')) ^ sin^ 



3 4 I . 2V' " ' V. 4 

2^ 



^2/1 „,2 I ^2 «^ 



-y [L-y +x 



\/A ^(1- cos 0)2 I (68) 



where 

A = A(l, 2/2, a;2) = 1 + + a;^ - 2y^x^ - 2y^ - 2x^ . (69) 
For the mf, 7^ normalized helicity fractions one now obtains 

Ql= {{l-x^ f -y\l + x^))/N , 

g+= y'il^y' +x^ -VX)/N, 

y\l-y^ +x^ + VX)/N, (70) 

where 

iV = (1 - 2-2)2 y2^^ _ 2y2 ^ ^2^ 

Let us compare the resulting numerical mi, ^ values for the normalized helicity fraction with 
their to;, = counterparts. One obtains (we take m\, = 4.8 GeV as default value) 

TOfc = 4.8 GeV: Ql-Q^- Q-= 0.7025 : 0.0004 : 0.2971 , 

TOb = : Ql--Q+-- Q-^ 0.7031 : : 0.2969 . (71) 
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The effect of including the nonvanishing bottom quark mass can be seen to be quite smah. 

Although wc have derived the decay distributions ((67|) and ([68| for the Born term case, the 
angular structure is quite general as will be shown in the next subsection. In the general case 
one has to replace the LO Born term structure fimction i/^''(Born) in (p7|) and by their 
generalized counterparts as e.g. the corresponding NLO or NNLO structure functions. 



3.2 Angular decay distribution for t b + W^{-^ + Ui) (II) 

The cos 9 dependence of L^^iJ^'' can also be worked out in a more systematic way by using 
the completeness relation for the polarization four- vectors Eq. One can then rewrite the 

contraction of the lepton and hadron tensors L^^H^'^ as 

= J2 L^^''''^^'{ni)el{m)el,{m')e,(m')H^^ 

ra.m' 



m.m 



We have thereby converted the invariant contraction Lf^^H'^'^ into a contraction over the spatial 
spherical components Lmm'Hmm' {m,m' = +,0,—), where the spatial spherical components of 
the lepton and hadron tensors are defined by 

H,nm' =Hf'''e*^im)e,{m'). (73) 

We have again dropped the g^g'' -terms in the completeness relation in Eq. since q^L^^ ^ = 
q^L^^ = for massless leptons. The nice feature of the representation ([7^ is that the left 
bracket and the right bracket in the next to last row of (|72p are separately Lorcntz invariant. 
One can therefore evaluate the left bracket in the rest frame, and the right bracket in the 
t-rest system without involving any boost. 

Let us now specify the the VF'''-rest frame four-vectors that are needed in the W^-rest frame 
evaluation of the lepton matrix L„im'- In the W'^ rest frame one has 

= miy/2(l; sin6'cos0, sinflsin^, cos0), 

= miy/2 (1; — sin 6* cos 0, — sin sine/) ,— cos 0) , (74) 

and the polarization vectors (in our convention = (oq, a) and = (oq, ~a)) 

e^{L) = (0; 0,0,-1), 

e^(±) = -i=(0;±l,z,0). (75) 
It is then straight-forward to evaluate Lmm' = L^'^ ef^{m)el(m') using the lepton tensor ([SO]). 



^ Since the method is general we can omit the LO specification in //^"(Born). 
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The various components of the lepton matrix Lmm' can be written in a very compact and 
suggestive way in terms of Wigner's small d^-function. One has 

L„„„. (9, 0) = m^^ dl 1 {9)dl, 1 (^)e'('"-™')^ , (76) 

where the spin one d^ function is given by (convention of Rose) 

/ i(l + cos6i) -^sin6' i(l-cos6i) \ 
dl^'iO)= ^sin0 cos9 -^sin0 . (77) 

y i(l-cos6i) ^sin6' i(l + cos6i) J 

The rows and columns arc labeled in the order (+1,0,-1). The representation ([75]) should 
be of no surprise to anyone who is familiar with the behaviour of angular momentum states 
under a rotation by the angles and (f>. In the lepton system {x',y',z') the only nonvanishing 
component of the lepton matrix is L^i^^i = as the antilepton and the neutrino are 

both left-handed (see Fig. 2). Eg. ([75]) represents the rotation of the lepton matrix from the 
lepton system (x' , y' , z') to the hadron system (x, y, z). In the case m; ^ one has to augment 
Ea.(|76[) by temporal spin components and interference contributions of the temporal spin 
and spatial spin 1 components |44| . 

When integrating L^i,H^'^ over the azimuthal angle (f> one remains only with the three 
diagonal elements of Hmm' ■ One has 

f dcj, L^.H^^" = 2nml. ^ <,+i(^)rfL,+i(^)^fmm 

m= + l,0,-l 

^2TT^m'^(^^sm^9Hoo + ^il + cosefH++ + ^{l-cos9)'^H^^^ . (78) 

By convention one drops one of the double indices in the diagonal elements of the hadronic 
density matrix Hmm, i-C one replaces Hqq — > Hq and H±± H± as has been done in the rest of 
this paper. For the LO case one reproduces Eq. (PT|) using Hqq{= Hl) ~ iqP + \H_i.iq\^, 

i/++(= il+) = and iJ__(= 7?_) = from 0. 

The advantage of method II is that the method can easily be applied to more complex 
decay processes involving spin. Also one can easily incorporate lepton mass effects and include 
polarization effects of initial and final state particles [33]. For example, method II was 
applied to the full angular analysis oi B ^ D , D* + i + vg = e, t) (331 [331 ^''^'i the rare 
decays B —f K, K* (^ = e, r) [45] including results on the polarisation of the final 

lepton. Another example is [42] where we have used method II to describe the semileptonic decay 
process of a polarized "EP , S°(t) Ti^ -i-vi {l~ ~ e^,fi^) followed by the nonlcptonic decay 
S+ p + tt'^. In this process the mass difference M30 — M-^+ = 125.46 MeV is comparable 
to the /x-mass which makes inclusion of lepton mass effects mandatory. In fact one finds 
r(/i)/r(e) « 1/120 in this process. A cascade type analysis as used in the method II is 
ideally suited for Monte Carlo event generators that describe complex cascade decays involving 
particles with spin. In fact, we wrote a Monte Carlo generator for the above semileptonic "EP 
decay process [42] which was profitably used in the analysis of the NA48 data on this process. 

3.3 Experimental results on helicity fractions 

An early MC study quotes experimental sensitivities of SQl = 0.7% and SQ+ = 0.3% for 
an integrated luminosity of 100 fb~^ at Tevatron II energies which corresponds to « 8 • 10^ 
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(ti)-pairs [36]. Compare this to the NLO QCD changes SQl = 0.7% and Sg+ = 0.1% to 
be discussed later on which shows that the radiative corrections are of the same order as the 
experimental sensitivities. Much higher event rates can be reached at the LHC in one year. 
A more recent MC study based on 10 fb~^ at the LHC quotes measurement uncertainties of 
SGl = 1.9%, dg+ = 0.22% and 60- = 1.8% [47]. 

Experimentally, there has been a continuing interest in the measurement of the helicity 
fractions. Latest measurements are 

CDF(2008 gH]) : Gl = 0.66 ± 0.16(stat) ± 0.05(syst) 

g+ = -0.03 ± 0.06(stat) ± 0.03(syst) 

DO(2009 gg) : Gl = 0.490 ± 0.106(stat) ± 0.085(syst) 

g+ = -0.104 ± 0.076(stat) ± 0.066(syst) (79) 

All of these measurements are well within the SM predictions. 



4 Construction of covariant helicity projectors 

In Ea. (j73p we have defined the helicity structure functions Hm [m = L, + , —) which multiply 
the angular factors in the angular decay distribution. According to their definition in Eq. (|75|) 
the helicity structure functions Hm can be calculated in a frame-dependent way by use of the 
frame-dependent polarization vectors (j29p . It is much more convenient to calculate the helicity 
structure functions covariantly, and, in fact, a covariant projection is indispensable for the 
NLO and NNLO calculations. The covariantization is achieved by defining covariant helicity 
projectors TP^ which covariantly project onto the helicity structure functions via 

This definition holds for any general hadron tensor structure irrespective of the fact that we 
have dealt only with the Born term hadron tensor up to now. To construct the covariant helicity 
projectors we start with their representation in terms of the t-rest frame polarization vectors 
according to the definition Eg. ([75]) . One has 

Pf ^ e*^{L)e'^{L), 

Pf = e*'^(±)e"(±). (81) 

In covariantizing the forms (|8ip it helps to remember that the helicity projectors must be 
four-transverse to the momentum of the , i.e. they must satisfy 

g^P^ = g.Pr = . (82) 

Further, they must satisfy the orthonormality and completeness relations 

orthonormality : g^i^PJ^j^ = — 1 

n fiW^^W'^" — - f5 P^" 

completeness : ^TP'^^ := P^,";^ = -y^" + ^ (83) 
As it turns out the covariant projectors can be constructed from the following three projectors 
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• Projector for the total rate Pj*/^^ 



K^L = -9''' + ^ (84) 



• Projector for the longitudinal helicity rate P^"^ 



• Projector for the forward-backward asymmetric helicity rate = — 
(eoi23 — 1) 

= —7^^e'""'^Pt.o.qf} (86) 
mtlql 

The denominator factor |(fp refers to the top quark rest frame. In invariant form the nor- 
malization factor is given by |(f p = {{ptqY' — mY/rnl) / ml . Finally, the three projectors read 
(iJ™ = Pj^r^^M- m = L,+,-) 



P 



pi/ _ '"-w 



^± = \ {^u+L - IPf ± IPr) ■ (87) 

It is instructive to check that, in the i-rest frame or in the W^"'"-rest frames, the covariant helicity 
projectors in Eq. (|87p reduce to the form ()8ip in terms of the rest frame polarization vectors 
([29|) and ([75|) . respectively. Note, though, that in the ly-rest frame the normalization factor 
|(f I in Eqs. ([55]) and ([55]) has to be replaced by y\pt\ where \pt \ is the magnitude of the top 
quark momentum in the iy+-rest frame. 

The denominator factors \q\~'^ and \q\~^ in Eqs. ([55)1 and (|86p are needed for the correct 
normalization of the projectors, cij. Ea. (j83[) . As we shall see later on the denominator factors 
and \q \ somewhat complicate the NLO and NNLO calculation of the helicity rates as 
compared to the total rate. 



5 Narrow width approximation 

Let us begin with by discussing how to factorize of the three-body rate V{t b + £^ + Ui) into 
the two-body rates r(i b+W'^) and r(iy"'" — > i'^ + i/i) using the narrow width approximation 
for the ly-boson. The rate formula for the three body decay t ^ b + + ly reads (see [27j ) 

which we write as 

Ta = ^i?3 [iMsH . (89) 
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The squared three-body matrix element jMsp is given by 



1 



W 



(90) 



where we have introduced the Breit-Wigner hne shape to account for the finite width of the 
H^+-boson. We have also reinstituted the factor of 8 • 8 = 64 in (|90p which was introduced 
earlier for convenience. 

Next we introduce the identity 



d\ 

2 



(91) 



which can be seen to be true in the rest frame where = E'^ and / ^g"! ^{Ew — Ei+ — 
E^i) = 1. The identity ([9T|) allows one to factorize the three-body phase space integral R^{t 
b + + vt) into the two-body phase space integrals i?2(i ^ & + W^) and R2{W^ vi). 
One has 



i?3 — 2mw 



df_ 

(27r) 



1 

2mt 



1 d^pb f 1 d^q 
(27r)3 2^ J '{2nfJE~ 



i2n)^S^^\pt-pt-q) 



2mw 



1 d^p, 



1 d^p^ 



(27r)3 2Eg+ J (27r)3 2E, 



■{2n)H^^\q-p,^-p,,))'j 



(92) 



The phase space nicely factorizes. But how about the factorization of the squared three-body 
matrix element jMap? The matrix element squared also factorizes after angular integration 
which can be seen by using the relation 



d cos 9 d4)Lmn {9, 0) 

which follows from the explicit representation of L, 

1 



47r 



(93) 



given in Eg. ([75]) . In fact, one has 



j d cos 9 d(j)(^^^ Hmn ^ Lr. 



Tr{(5„i„} 



d cos 9 ( 



(94) 



The factor 1/Tr{(5mn} = 1/3 provides for the crucial statistical factor \/{2sw + 1) in the W 
width formula. Note that the explicit angular integrations over cos 9 and cf) appearing in 
are implicit in ((HH) Q- One thus finds 



3 



dq^^R2 \M\^{t^b + W^ 



\M\^{W+ ^ e+ + iy) 



1 



q^ -mi^+imw^W 



(95) 



""We mention that an alternative derivation of the appearance of the statistical factor 1/3 has been given in 
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The narrow-width approximation consists in the replacement of the Breit-Wigner hnc shape 
by a (5-function, cif. 



1 



-5{q^ - m^) . (96) 



Using the narrow- width approximation for the iy+-boson the three-body decay t ^ b + i'^ + v 
can be seen to factorize, cif. 

Tw 

= T{t^b + W+) BR{W+ + vi) (97) 

which is a result which one expects from physical intuition. Incidentally, the derivation of 
the factorization formula (|97|) was posed as one of the problems in the 2004 TASI lectures 
of T. Han |51| . Judging from the contents of this subsection this was not one of his simpler 
problems. 

The numerical value of the finite-width correction to the total width listed in ([6]) consists 
of the replacement of 5{q^ — m^) by the Brcit-Wigncr line shape and integrating over , cif . 



dq'd{q'-rnl.) ^ dq'^^^^- , , , (98) 



(g2 - m^^.f + rn^^Tl 

where Tw is the width of the W-boson {Tw =2.12 GeV ). 

Numerically the finite- width correction to the total Born term rate amounts to 1.56% (see 
Eq.®) and is of the order of Tw/'mw — 2.64% as would be expected. A more extensive 
discussion on finite width corrections can be found in |50[ 152] . 



6 Higher order corrections to helicity fractions 
6.1 NLO QCD and electroweak corrections 

As in the calculation of the NLO total rate structure function ffy+L = + i7_ + Hl (we 
call -\- H- = Hjj where U stands for the "unpolarized transverse" ) we have employed the 
traditional technique when calculating the helicity structure functions and H±^ i.e. we have 
separately calculated the hadronic loop and tree contributions after contracting them with the 
relevant projectors P^'', . 

As mentioned before, the appearance of the normalization factors \(f\^^ and |(f|^^ in the 
projectors make the calculation technically more difficult than that for the total rate. For the 
one-loop contribution the additional normalization factors are of no concern since they appear 
only as overall factors outside of the one-loop integral. This is different for the phase space 
integration of the tree-graph contributions where the normalization factors appear under the 
integral. Typically one of the phase space integrations is over the scaled invariant mass of the 
bottom quark and the gluon z = (pf, + pg)'^ /m^. The normalization factors then appear as 
overall factors \q\~^ and \q\'''^ in the phase space integral, where 

k1 = ^VA(l,y2,^). (99) 
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The ensuing class of phase space integrals is more general and more difficult than the class of 
integrals appearing in the total rate calculation. Nevertheless, the phase space integrations can 
still be done in closed form. 

As a sample result we present the mi, = result for Tl = ri/r[/+L(Born). One obtains [HI 

MM 

- ^(l + 52/2 + 2/)-3(l-y2)2 ln(l-y2) + 16i/{l + 2y^)\n{y) - 2{l-yf 

X (2^y + 62/2 + y3) In(l-y) \n{y) - 2{l + yf{2 + y + 6y^'-y'') \n{y) ln(l+y) 

- 2(l-y)2(4+3y + 8y2+y3)Li2(y)-2(l+2/)2(4-3y + 8y2_y3)Li2(_y)|) . (lOO) 

In the limit y ^ one finds f l(NLO) = 1 + a,CF/{2n){b/2 - 2tt^/3) and thus f l(NLO) 
saturates the total rate r(NLO) (see Ea. ([57|) ) in this limit. This is expected since Tjj/Tl oc 
ra^/ml = . Results for the other two NLO QCD helicity rates r+ and r_ can be found in 
[HI mi I3S] ■ The NLO electroweak corrections to the helicity rates can be found in [3] . 

Let us summarize our numerical NLO results on the helicity fractions including also the 
finite width corrections discussed in Sec. 5. We write 

r, = r,(Born) + Ar,(QCD) + Ar,;(EW) + Ar,(FW) + Ar,(TO6 ^ 0) . (101) 

As before we normalize the partial rates to the total Born term rate r[/+i(Born). Thus we 
write Ti = ri/r[/+L(Born) (i = +,— ,L). For the transverse-minus and longitudinal rates we 
factor out the normalized partial Born rates F^ and write (i = — , i) 



F, = F,(Born) [l + (5.(QCD) + 5,{WN) + (5.(FW) + 5T,{mb ^ 0)j , (102) 

where 5i = r^z+i (Born) AF^/Fi (Born). Writing the result in this way helps to quickly assess 
the percentage changes brought about by the various corrections. 
Numerically one has 

f_ = 0.297(1 - 0.0656(QCD) + 0.0206(EW)-0.0197(FW)- 0.00172(mb 7^0)) 

= 0.297(1 - 0.0664) , (103) 

and 

Tl = 0.703(^1 - 0.0951(QCD) + 0.0132(EW)-0.0138(FW)- 0.00357(TOfa 7^0)) 

^ 0.703(1 - 0.0993). (104) 

It is quite remarkable that the electroweak corrections almost cancel the finite width corrections 
in both cases. 

In the case of the transverse-plus rate the partial Born term rate cannot be factored out 
because of the fact that F+(Born) is zero. In this case we present our numerical result in the 
form 

f + = Af +(QCD) -f Af +(EW) + Af +(mb 7^ 0). (105) 
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One has 



f + = 0.000927(QCD) + 0.0000745(EW) + 0.000358(mfc ^ 0) 

= 0.00136. (106) 

Note that the finite width correction to the transverse-plus hehcity rate is zero. Numerically 
the NLO corrections to r+ occur only at the pro mille level. It is save to say that, if top quark 
decays reveal a violation of the SM left-chiral iy — A) current structure that exceeds the 1% 
level, the violations must have a non-SM origin such as e.g. an admixture of a right-chiral 
{y + A) current structure in the decay vertex t ^ 6 + 



6.2 Quality of the (rnvtz/m^)— expansion 

In order to check on the quality of the y = (TOn^/mt)-expansion we take the known closed form 
NLO result (|100p for and expand it in powers of and y^lny. The expansion of the curly 
bracket in (jlOOp reads 




(107) 



Note that TL{as) as the phase space closes at y = 1 (Li2(— 1) = — •7r^/12). In Fig. 7 we 
show a plot of the y-dependence of Fi,(as) (in units of [Fq ^Cf]) for diS'erent orders of 
and for the full result. All curves start at (5/2 - 2tt^/3) = 0.459 for y = 0. The full result 
goes to zero at y ~ 1 remembering that Li2(— 1) = — 7r^/12. As Fig. 7 shows the quality of the 
expansion is already quite good at 0{y^) even for large y- values. 

This raises the hope that such a (mvK/'T^t)-expansion can also be usefully employed in other 
contexts. One could think of possible applications of the NNLO calculation oi t ^ b + VF+ 
discussed earlier (which only exists in expanded form) to processes such as 

• b ^ u + £^ + vi 

• ^ + + l^fi 

extending over the whole kinematical range 0^5^^ (mi — m2)^ in these processes. 

The region very close to the upper kinematical limit of given by q^^^ = {nit ~ mi,)'^ 
requires a separate discussion because this region is sensitive to mf, 7^ effects. The upper 
kinematical limit is called the zero recoil point since g = at this point. For example, at zero 
recoil (y = 1 — mh/mt = 1 — 4.8/175 = 0.973) one finds 

^ : QL:G+:Q- = l/i: 1/3 : 1/3 (108) 

using Eg. ([70)1 . The equipartitioned helicity fractions result from the fact that, close to zero 
recoil, the only surviving transition is the allowed Gamow- Teller s-wave transition. However, 
for mf, = one has the zero recoil ratios at y = 1 (see Eq.®) 

mh = : : e+ : = 1/3 : : 2/3 . (109) 
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0.2 0.4 0.6 0.8 1 

y=m^/m, 



Figure 11: Quality of the y-expansion of the as correctious to TL{as). Shown are different 
orders of the y-expansion in units of [Tq as Cf/2tt]. Dashed hne: 0{y'^); dotted hne: 0{y^); 
dashed hne : O(y^); dash-dotted hne: 0{y^); full line: exact result. Vertical line corresponds 
to the physical point y = mw/rnt = 0.459 . 



In order to investigate the behaviour of the helicity fractions close to zero recoil, in Fig. 8 we 
plot the y^-dependence of the helicity fractions for mi, ^ and rrib = with zero recoil values 
at y = 1 — rrib/mt and y = I, respectively. In the region close to their respective zero recoil 
points the curves considerably differ from each other. Away from zero recoil the rrib = and 
TTib 7^ curves very quickly approach each other. Fig. 8 shows that it is safe to use the = 
approximation for y-valucs below y w 0.9. 




Figure 12: Helicity fractions Gl,G+ and G- close to zero recoil. Dashed line: nib 7^ 0; full line: 
TTib = 0. The kinematical zero recoil point is given by y = 1 (full line) and y — 1 — nib/mt 
(dashed line). 
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6.3 NNLO QCD corrections to helicity fractions 



In Sec. 2.4 we have desribed how the total NNLO rate can be calculated in a y = {mw/^t)- 
expansion using the optical theorem. Two new features appear in the corresponding NNLO 
calculation of the helicity rates Tl,±- First there is a parity violating three- loop contribution 
which is projected out by the projector Wp^. One has to deal with the problem of how to treat 
75 in the environment of dimensionally regularized loop integrals. We take the prescription 
of [S3] and replace 

7m75 ^ ^W77"7^^ • (110) 

When using this prescription one needs to add finite three-loop counter terms which arc given 
in Hi]. 

The second new feature is related to the normalization factors \q\^^ and |(f in the three 
helicity projectors P^'^j. which replace the total rate projector IP[5^^ = — g'^" -I- q^q'^/m^. 
In the hard region one can expand in inverse powers of the (large) propagator pole factor 

P={Pt + q? - ml 

2 

-m^^. (Ill) 

One expands in the propagator pole factor P = {pt + q)"^ — mf = 2ptq + q^, i.e. ptq = 
ip(l — m'^P~^) where one can replace q^ by to^ since one is cutting through the M^-line 
anyhow. One then has 




p2 



n=Q 
oo 



2m^F2 _ 



1 2mt \ ^ 



n=0 



p2 



4P2 



(112) 



Thus, the additional propagator-likc structures from the projectors arc transformed into a scalar 
on-shell propagator with momentum p + q and mass mt raised to arbitrary, integer powers. This 
will eventually lead to twelve additional three-loop master integrals next to the master integrals 
appearing in the total rate calculation of [TH [20] whose imaginary parts can again be calculated 
in closed analytical form using the cutting rules. 

In the soft region one cannot perform an expansion of \q\, since \q\'^ ~ Qo ^ ""^"w 90 is of 
order mw in the soft region. However, in this region the W boson loop factorizes. Therefore, 
one only has to replace the usual one-loop vacuum bubble integrals with integrals of the type 

Iif-mi,)ili~m^^r ' ^^^^^ 

with n = 1/2 and 1. These integrals are not difficult to evaluate. 

The validity of the treatment of these two new features has been tested against the known 
NLO results up to 0{{mw /iritY^) [35]. First results of the NNLO calculation have been 
published in [38]. Complete results on the NNLO calculation of the helicity rates will be 
published soon [55] . 
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7 Summary and conclusions 



We have discussed some of the properties of the top quark with an emphasis on the SM decay 
properties of the top quark. We have defined partial hehcity rates into polarized M^"^-bosons 
and have derived the resulting angular decay distribution of W^^^ l^v. We have described 
the LO calculation of the partial helicity rates using several methods including also the optical 
theorem and a mw /mt-expansion as a preparation for the description of the NNLO calculation 
of the total rate and the partial helicity rates. We have summarily described the main features 
of NLO QCD and electroweak corrections to the total width and the partial helicity rates. 

We are looking forward to the LHC era with its expected wealth of data on the top quark 
and its decay properties. 
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